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Abstract 

We consider ring polymers in good solvents in the dilute limit. We determine 
the structure factor and the monomer-monomer distribution function. We com- 
pute accurately the asymptotic behavior of these functions for small and large 
momenta and distances by using field-theoretical methods. Phenomenological 
expressions with the correct asymptotic behaviors are also given. 



1 Introduction 



The statistical properties of dilute polymers in good solvents have been the subject of 
extensive studies during the years A significant understanding of the problem 

was reached when it was realized that long polymers could be modelled by chains 
with an excluded-volume interaction. This allowed the introduction of simplified 
theoretical models which could be analyzed more easily. From a theoretical point 
of view, an important step forward was made by de Gennes 0, who proved that 
the statistical properties of dilute polymers could be obtained as the limit n ^ of 
the n-component (j)'^ theory, opening the field to the many methods that have been 
developed for the study of the critical behavior of spin systems. 

In nature polymers may have many different geometrical conformations. In this 
paper we will focus on ring polymers and we will determine the structure factor 
and the monomer-monomer distribution function. For this purpose we will use the 
mapping with the theory p-pHf and some recent theoretical results for the two- 



point function of biquadratic operators [11, 12|. Indeed, we show that the structure 



factor is directly related to the energy-energy correlation function. Then, by extending 
the results available for the end-to-end distribution function ||T3|- pO|] , we compute its 
large- and small-momentum behavior. Moreover, by using an interpolation formula 
based on a dispersive approach [|12|,|21[], we obtain a general expression valid for 



all momenta, i.e. for all values of the radiation (neutron) wavelength. We also 
investigate the monomer-monomer distribution function that is closely related to the 
structure factor by a Fourier transform. Its properties for small and large distances 
are investigated in detail. 

The paper is organized as follows. In Sec. |^ we define the structure factor and 
the monomer-monomer distribution function for A^-step ring polymers. In Sec. p.l| 
we show that the computation of the structure factor is equivalent to the calculation 
of the energy-energy correlation function in the n-component theory in the limit 
n 0. Field-theoretical results for such a correlation function are reviewed in Sec. 
3^ . In Sec. § we study the structure factor. In particular, we determine its small-g^ 
behavior extending the classical Guinier formula, its large-g^ behavior, and give a 
general interpolation formula that has the correct asymptotic behavior for — * 
and — > cxD. Finally, in Sec. ^ we report the corresponding expressions for the 
monomer- monomer distribution function. We also verify that the phenomenological 
expression often used for the end-to-end distribution function for linear polymers 



15-18 also provides a good approximation to the monomer-monomer distribution 



function considered here. 



2 Definitions 

We consider a monodisperse ensemble of ring polymers with N monomers, labelling 
the monomers from 1 to A^. Then, we consider two monomers i and j of the walk and 
the unnormalized distribution dij^jsilf) of the distance between the two monomers. 
Such a distribution function depends only on the relative position \i — j\ of the 



2 



monomers along the polymer. Then, we define its average over all pairs i.e. 



1 ^ 



Pn{^ = 7?I E ^^.,^(^1' (2-1) 
the corresponding normalized distribution 

P.irl = . (2.2) 

and the mean squared radius of gyration 

Rl^ = lJd''s\s\'Pr,{s). (2.3) 

On a lattice, dij^Ni'f^) can be identified with the number of iV-step lattice rooted self- 
avoiding polygons such that monomer i is in the origin and monomer j is in r. Note 
that 

Pn = Y1 dij,N{s) = ^^PNis) (2.4) 

s s 

is the total number of A^-step lattice rooted polygons. 
The structure factor is defined by 

SNiq) = I d'^se^^-'P^i^, (2.5) 

which, by definition, satisfies Sn{0) = 1. For elastic scattering, the momentum q 
is directly related to the wavelength A of the incoming radiation in the scattering 
medium and to the scattering angle 9 by 

, , 47r , , 

kl = ysm-. (2.6) 

For comparison, we will also consider linear polymers and correspondingly we 
define the mean squared end-to-end distance -ReAf the mean squared radius of 
gyration R^^n- In the limit ^ oo, the radius of gyration and the end-to-end 
distance diverge with the same critical exponent z/, i.e. 





= %N^-, 












= d,N^\ 


(2.7) 



In two dimensions the universal exponent v is given by v = 3/4, while in three 
dimensions 

^ r 0.5877(6) Ref. [22] 

10.58758(7) Ref. H. ^ ^ 



See Ref. \UM for other numerical and theoretical estimates. 
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In the limit N ^ oo, the normahzed distribution PNir) and the structure factor 



Sn^q) obey general scaling laws. More precisely, for iV — oo, \r\ 
p = f/ Rg^N and Q = qRg^N fixed, we have 

1 



oo, \q\ 



0, with 



SN{q) 



siQ). 



/(P), 



These functions satisfy the following normalization conditions: 

'dpf{p) = l Jdpp'f{p) = 2 

s{0) = 1 



dsiQ) 



(2.9) 



(2.10) 



3 Field-theoretical results 



3.1 Generalities 

We wish now to derive some general properties of the functions we have introduced 
above. We will use the Laplace-de Gennes transform method [B|,|0,|18i and the results 



for biquadratic correlation functions in n- vector models obtained in Refs. |lTT|,pr2 



For this purpose, consider a hypercubic lattice in d dimensions and the n-vector 
model Hamiltonian 

H = -[3Y.^r-(Ts, (3.1) 

(rs) 

where cr^ is a unit n-dimensional vector and the sum is extended over all lattice 
nearest-neighbor sites (rs). Then, consider the correlation function 



Cm(^; P) = ((cTo ■ cr^)(crr ■ cr,.+^)) - ((ctq ■ cr^))^ 



(3.2) 



where p is an arbitrary lattice unit vector. By using standard results (see, e.g., 
Refs. 0,H]), it is easy to verify that 



lim 

n-*0 



n 



N-2. 



(3.3) 



N 



where CN,^i{r) is the number of A^-step lattice self-avoiding polygons that go through 
the lattice links (0, p) and (r, r + p). In the limit — * oo, |r| — oo we are interested 
in, C]s[^^{f) is proportional to the number of lattice polygons PNif) going through 
and r. Then 

-C^(f; 13)] ^ const X P^^p^P^^r). (3.4) 

Now, let us consider the left-hand side. In the critical limit it can be identified 
with the energy-energy correlation function in the theory. More precisely, if 
= {Pc — P) / Pc ^ is the reduced temperature in the n-vector model, we have 



lim 



C^(f; p) ^ const X ^^(f, kt„), 



(3.5) 



4 



where is a constant and 

G^(f,t)^(02(O)0^(f)),-(02(O))?, 



(3.6) 



where {■)t indicates that the average should be taken at reduced temperature t in the 
theory. 

For iV — > oo one can then derive (see, e.g., Ref. M) 



PN^r) ^ const X 



— re hm 

ioo 2TCi n^O 



n 



-GE{r,kt) 



(3.7) 



where the integral is over a line parallel to the imaginary axis that leaves all singu- 
larities in the left-hand side. 



3.2 Behavior of Gsif, t) 
3.2.1 General results 

The general properties of GE{r,t) and of its Fourier transform GE{q,t) were studied 
in Refs. []TT1, |T^] . We will review them here. In the scaling limit t ^ 0^ , we have 



Gz,(g,t) = A+r°/^(|g|0 i + oir,t^' 



(3.8) 



where a = 2 — di> {a = 0.23726(21) in three dimensions), /^(y) is a universal function, 
and ^ is the second-moment correlation length defined in terms of the two-point 
function of the fundamental field G^{q, t), 



-G^(0,t) 



-1 dG^{q,t) 



(3.9) 



52=0 



Notice the presence of two types of corrections in Eq. ( |3.8| ). First of all, there are 
the corrections due to the leading irrelevant operator with exponent A. Additionally, 
there are corrections with exponent a due to the presence of a background term. Note 
that both in two and three dimensions p5[ a < A, and thus the leading correction is 
due to the background term. 

The function fE{y) has a regular expansion for y — > 0, 



2n 



(3.10) 



n=l 



while for ?/ ^ 00 it has the nonanalytic behavior 



fE{y)^E^y'^'- l+E^y 



(3.11) 



For generic values of the universal function fE{y) can be written in the form 



fE{y) = 1 



y'^Ei . f 7ra\ 



sin -— 



71 



\2v ) hs 



dx ■ 



X + y"^ 



-Fe(x) 



(3.12) 
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where Ei is the coefficient appearing in Eq. ( p. Ill) , is a constant, and Fe{ii) is the 
spectral function. The constant S~l.i is universal and it is defined by S'^ = ^^/^gap, 
where ^gap is the exponential correlation length that determines the large-distance 
behavior of the two-point function G^(r, t) of the field 0. The constant Ei and the 
spectral function are related by the sum rule 



El f-Ka 
— sm — 

TT \2v ) J 48 



poo 

/ c/xx-^-"/^^'^) FEix) = 1. 



(3.13) 



The representation (|3.12| ) is exact under rather mild assumptions [|12],^. Approx- 
imate expressions can be obtained by choosing approximate forms for the spectral 



function. As in Ref. 121, we choose 



(3.14) 



where E2 and E3 are the constants that parametrize the large-g^ behavior in Eq. 
( CT ) and 



$2 

$3 



7r(l — a) 7r(l — a) ira 
cos h sm cot — , 

71 71 Tca 

COS h sm — cot — . 

2u 2iy 2iy 



(3.15) 



We further define 
g{p; x) = 



TT 



sm 1 1^ + Tip 



X 



-l-p-a/{2u) 



+ 



2v (45+ ) 



l,l;l-p- — 



a 45 



{a^2pv){x^^SXi) V 
where 2-^1 (a, c; z) is a hypergeometric function [^]. Then 



9 ' '(3-16) 
2v x + ASlJ 



fEiy) = 1 - 



■ sm , 
TT \2iyJ 



g{0;y^)+E2<!>2g{K-^-^y^ 



or, by using the sum rule ( |3.13| ), 



fE{y) = — sin (^) 
7C \2v J 



2v 



\ — OL 



E;^,9\-,y' 



(3.17) 



2v 



2v' 



From Eq. ( p.8| ) we also obtain 



OL jA> 



— fE{s) l + 0(t",t 



where s = r/^, and 



{2tiY 
6 



/E(y)e^^X 



(3.18) 
(3.19) 

(3.20) 



Table 1: Estimates of the small- momentum expansion coefficients by using the fixed- 
dimension expansion in d = 3 {d = 3) and the e expansion (e-exp). "final" labels our 
subjective final estimates. Perturbative series from Ref. |]12|. 





ei 


62 


63 


64 


65 


id = 3) 
(e-exp) 
final 


-0.0322(11) 

-0.0323(7) 

-0.0323(10) 


0.00382(15) 
0.00398(11) 
0.00390(15) 


-0.000597(23) 
-0.000636(10) 
-0.000620(30) 


1.06(4)xl0-^ 
1.14(4) X 10-"^ 
1.10(8)xl0-'^ 


-2.03(8) X 10-5 
-2.22(3)xl0-5 
-2.15(15)xl0-5 



For s — >• 00, fsis) decays exponentially as 



where 



2\/ o « 




(3.21) 



(3.22) 



and p is an exponent that can in principle be computed perturbatively. 
An approximate expression can be obtained from Eq. (|3.18|) : 



fE{s) 

where 



El . /Tca 
— — sm — 



1 — n 

h{Q-s)+E2^2h[^—;s 



2u' 



h{p; s) = 2s2f-2+°/'^r [2-2p- a/u, 2s^Sij i , 
and r{a;x) is an incomplete F-function p6 |. 



(3.23) 
(3.24) 



3.2.2 Numerical results 

We have introduced in the preceding section several constants that will be computed 
here. 

First, let us consider the constants e„ that parametrize the small- momentum ex- 
pansion of fEiy)- Perturbative series were derived in Ref. []T2[ for e^, n = 1, . . . , 5, 
in the fixed-dimension expansion in three dimensions (four loops) and in e-expansion 
(three loops). We resummed the perturbative series by using their large-order behav- 
ior and performing a conformal mapping as in Ref. Mean values ad errors were 
computed using the algorithm of Ref. . For the expansion in fixed dimension we 
used for the four-point renormalized coupling |]29| g* = 1.395(15) that includes all 
available estimates [p0|-[3^. The results are reported in Table |l[ Note that the two 
expansions give estimates that agree within error bars, confirming the correctness of 
our results within the quoted errors. 
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We can also use the results of Ref. |0 to compute the large-momentum expansion 
coefficients Ei. In e = 4 — d expansion, they are explicitly given by 



El = l + l + Oie'), 
E, = -2 + | + 0(a, 

Es = 2-^ + 0(e2). (3.25) 

In the following we will also be interested in the product E1E2 given by 

EiE2 = -2 + ^ + 0{e^). (3.26) 

The large size of the coefficients makes it difficult to resum the perturbative series. 
We report here the results obtained by setting e = 1 and as error we quote the last 
coefficient: Ei = 1.5(5), E2 = -0.75(1.25), E^ = 0.25(1.75), E1E2 = -1.75(25). For 
the sum E2 + -E3 an additional term is known: 

E.^E.^-l.^(f^-'-).' + Oi.\ (3.27) 

SO that E2 + E^ = -0.49(1). 

We also need the universal ratio S'^f. In three dimensions it has been estimated by 
a variety of field-theoretical and exact-enumeration methods [Q, obtaining S*^^ — 1 = 
—3(1) X 10~^. In two dimensions, an exact-enumeration study gives ||3^ Sj^.^ — 1 = 
1(2) X 10"^ 



We wish finally to determine the functions fEijj) and /£;(p) by using Eqs. ( p.l8|) 



and ( |3.23| ). For this purpose, we must fix a, z/, S'^^, i?25 and E^, while Ei is fixed 
by using the sum rule ( |3.13| ). For v and a = 2 — 3z/, we use p3[ v = 0.58758, 



while for S]^^ we use the above-reported result. For E2 and E^, the only available 
estimates are those obtained in the e-expansion approach. They have a large error, 
while, apparently, their sum is more precisely determined. We have thus followed the 
following strategy. We have fixed E2 + E^ = —0.49 using the e-expansion. Then, we 
have chosen E2 and E^ so that for small y the approximation gives 1 — 0.0323?/^, in 
agreement with the results of Table |l[ In this way we obtain, E2 = —1.38, E3 = 0.89, 
and El = 1.60. These estimates are in good agreement with the e-expansion results, 
confirming that our approximation is reasonably correct for y 00. Moreover, it also 
nicely reproduces the small ^/-behavior. Indeed, it gives 62 ~ 0.0040, 63 ~ —0.00065, 
64 ~ 1.2 X 10~^, which are close to the results of Table |I|. 



4 The structure factor Sn^q) 

Using Eqs. (|2^), (p. 71), and ( |3^ ) we can determine the structure factor Sn{(i) in the 
scaling limit |g| — > 0, — 00, with Q = \q\R]y^g fixed. Neglecting corrections of 



8 



order N see the discussion in Sec. 3.2.1, we can write — from now on, the hmit 



n — > is always understood — 

/+i00 rjf . 



(4.1) 



where we have written for the second-moment correlation length ^, ^ ~ x^t''^ for 
t ^ 0. The nonuniversal factor can be eliminated by introducing the end-to-end 
distance for linear polymers. A simple calculation gives RIlB] 



where 7 is a universal critical exponent. If we define 

1 r(7 + 2v) 



(4.2) 



then 



2d r(7) 

dt 



(4.3) 



-ioo 2m ^ ' 

We can now use the results of the preceding section. First, for |g| ^ we obtain 



S'iv(g) = 1 + ^ e„ 



\ V{a + 2nv) 



(4.5) 



In order to express Sn{(i) in the scaling form (|2.9|) , we must express R^ j^ in terms of 
R^^N- Using Eqs. ( ^.31 ) and ( |4.5| ) we have in d dimensions 



2 _ , dSN{q) 
Rg^N- d 



-dei „, 1 , kR^ 



(j2=0 



T(a + 2u) 



(4.6) 



which allows to determine the ratio 

H 

Setting 



^--de.n^ 

^ — — Ucl r\j 



T{a + 2v) 
1 r(a + 2z/) 



dei r(a) 



we have 



5'Ar(g) = 1 + ^ e 



Via) 



\ r(a + 2nv) 



(4.7) 

(4.8) 
(4.9) 



In three dimension {d = 3), using ^ u = 0.58758(7), a = 0.23726(21), and |6 
7 = 1.1575(6), we obtain 



K = 0.21315(12), 



R = 2.40(8). 



(4.10) 
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Thus 



and 



= 0.089(3), (4.11) 



SN{q) = 1- -g2 + 0.061(5) X g^- 0.0073(7) X 
3 

+0.00063(10) X g^ - 0.000045(8) x Q^^ + 0{Q^^) 



(4.12) 



Note that the coefficients decrease quite rapidly, so that this expansion provides a 
good approximation up to g ~ 2. 

We can compare the resuh (|4.11|) with the existing estimates. The amphtude Sg, 
cf. Eq. ( |2.7| ), is quite well determined for self-avoiding walks on the cubic lattice: 
Ref. quotes de = 1.2167(50), while Ref. gives ~ 1.225. The amplitude ttg, 
cf. Eq. ( ^.71 ), for self-avoiding polygons on the cubic lattice can be obtained from 
the data reported in Ref. Fitting their data for R^ p^ with agN"^" + bgN^"'^ and 
using the estimates of u and A of Ref. we obtain ~ 0.102. Unfortunately, no 
error bars are reported in Ref. |3^ and thus we cannot quote an error on our result. 
Moreover, the algorithm of Ref. |^ samples only polygons with trivial knot type, 
and thus some systematic error can in principle be present. Using these estimates, 
we obtain H ^ 0.083, which is reasonably close to our result ( [4.11|) . 

It is interesting to compute ei and R in two dimensions, by using the available 
estimates oi H, H = 0.14605(7) (Ref. ^) and H = 0.1459(2) (Ref. O). Since 
u = 3/4 and g 7 = 43/32, we obtain d = -0.08447(4), R = 3.3398(16). 

For g^ ^ cxD, we can use Eqs. ( [4.1| ) and ( p.ll| ) to obtain 



SNiq) ~ (a - 1)^1^2 {f^q^RlN 



-l/{2u) 



(a - l)EiE2R-^^^^'^Q-^^''. (4.13) 

In three dimensions, using Eq. ( |4.10| ) and the estimate of E1E2 of Sec. |3.2.2| , we have 

SN{q) ~ -0.363(10) x EiE2Q'^^'' ^ 0.63(9) x Q-^^\ (4.14) 

Finally, we can determine an approximate form for Si\i{q), by using the interpolation 
formulae ( ^.17[ ), ( |3.18| ) derived in Sec. |3.2.1| . We obtain 



SN{q) 



1 + (a 
El 



sm 



71 



l)EiE2R-^'^''Q'^''' 

mQ^) + E2^2g 



(4.15) 



2^ 



where 



1v 



-p-al{2v) 



a + 2pz/ 

X2i^i ( 



-Via) 



1 — a 
1v 



dt 



■Q']+Es<!>3g(—;Q 



RQ'^t 



-a-2u 



a 



X 



2z/' 4Sli + RQH 



(4.16) 



The evaluation of Siy{q) requires the evaluation of the integral ( [4.161) . Particular 
care should be given to the branch cuts. There are indeed three cuts in the complex 
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a 0.6 



0.4 



0.2 





large q 




small q 




Guinier 







0.1 



10 



100 



Figure 1: Structure factor Sn{(1) versus = q^Rg. We also report the large- and 
small-Q^ behavior, cf. Eqs. ( [4.13D and ( [4.12[ ), and the classical Guinier formula 
SN{q} = 1 - QV3. 



t-plane: the negative real axis and the lines t = aexp[±i'7r/(2z/)], where a is real 
satisfying a > [rQ"^ / {4:S^)Y^^'^'^\ In practice, we have found convenient to integrate 
over the lines t = aexp[±z7r/(2z/)] + b, where a runs over the positive real axis and b 
is a fixed positive constant that we have taken equal to one. 

Using the values of the parameters determined in Sec. j.2.2 , we obtain the curve 
reported in Fig. |l|. We also report the small-Q^ expansion (|4.12| ) and the large- 



Q expansion ( [4.13D . Note that the small-Q approximation is indistinguishible in 
the graph from the full curve up to ~ 5, while the large-Q^ approximation is 
reasonably accurate for > 40. Numerical estimates for intermediate values are 
reported in Table 0. 

We can compare our field-theoretical determination with the numerical results of 
Ref. [0. In Fig. || we plot Sn{(i) in such a way that it can be directly compared 
with their Fig. 14. We observe a nice quantitative agreement. 



5 Monomer- monomer distribution function Piv(r) 

We wish now to compute the distribution function P/v(r), or, more precisely, the 
scaling function /(p) defined in Eq. ( |2.9| ). For this purpose we can either use Eq. 
( ^75|) and the results of Sec. § for Sn{(i), or use directly the field-theoretical results 
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Table 2: Estimates of the structor factor S'jv(g) for 4 < < 40, obtained by usinj 
the interpolation formula ( |4.15| ). 





SN{q) 




SN{q) 


4 


0.305 


16 


0.0733 


5 


0.241 


18 


0.0667 


6 


0.196 


20 


0.0612 


7 


0.164 


24 


0.0527 


8 


0.141 


28 


0.0464 


9 


0.124 


30 


0.0437 


10 


0.111 


34 


0.0393 


12 


0.093 


38 


0.0358 


14 


0.082 


40 


0.0343 
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available for the scaling function fE^p), cf. Eq. (|3.20|) , and 

rioo rjf ~ / , s 

/(p) = V{a)-K~'^'^ £^ 1^ eH-'^^^'^ fE [pt^K'l') . (5.1) 
Using the small-Q^ result ( [4 .91) , we can compute the moments of /(p), i.e. 

M, = I dVlpl'V(p). (5.2) 

A simple computation gives 

M, = {-2Kfk\e. \l n + 2n). (5.3) 

Numerically, in three dimensions, we have M\ = 2, M2 = 7.3(5), M3 = 37(4), M4 = 
230(35), M5 = 1780(330). 

We can use the large-g^ expansion to compute the small-p behavior of /(p). Ex- 
plicitly, by using Eq. ( [4.13|) , we obtain 



/(p)^po//^-^ (5.4) 

where 

Po = ia- l)E,E,(4^)-V(^-).-^/^ ^(^/^-^^/(^'-)) . (5.5) 
Numerically in three dimensions, Pq = —0.0250(7) x E1E2. If we use the estimate 



E1E2 = -1.75(25) reported in Sec. |3X|we have po = 0.044(6). 

Finally, we determine the large-p behavior of /(p). Starting from Eq. ( |5.1| ) and 
using Eq. ( p.21|) , we have 



f{p) ^ V{a)A-K-'"'-^/'p^ _e*t-+-'^+-Pexp ( -2pr^SlJ-K ] . (5.6) 



The calculation is analogous to that presented in Refs. [0,|T8[. We obtain that /(p) 
obeys the Fisher law 

/(p)~p'^exp(-Dp^), (5.7) 

where 

6 = (5.8) 

D ^ lAl^J^f, (5.9) 



and q is an exponent that can be expressed in terms of p, cf. Eq. ( 3.21 ). Numerically, 
in three dimensions, 5 = 2.425(4), D = 0.360(14), while in two dimensions 5 = 4, 
D = 0.15132(16). 
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large p 
small p 
phenomenological 




1e-05 



Figure 3: Monomer- monomer scaling function f{p) versus p. We also report the small- 
p behavior (|5.4| ), the curve 0.036e~^^ (labelled "large p") that has the correct asymp- 
totic behavior for p oo, and the phenomelogical representation 0.060p^^'^~^e~^'^'^^P . 



Finally, we can use the interpolation formula ( p.23| ) to obtain f{p) for generic 
values of p. We write 



El ( T\a 
— -— sm — 
47r2p \2v 



/i(0; p) + E2^2h ( + E^^^h p 



(5.10) 



2r(a) / 2 ^\p-l+a/{2v) r+i°° dt 



fip) 

where 

h{p; p) 

(5.11) 

Using the parameter values determined in Sec. |3.2.2| we obtain the curve reported 
in Fig. |[ 

In three dimensions the end-to-end distribution function for linear polymers is 



{p'K- 
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well described by a very simple phenomenological expression [p!5|-p!9 

/ph(p) = Aphp''e-^^''^\ 



(5.12) 



where 5 and 6 are fixed by the large-p and small-p behavior, and the constants A^^i 
and Dph by the normalization conditions ( |2.1CI| ): 

r[(i-z.)(^^ + 5)]^'/^ 



^ph 
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5_ 

47rr[(l-z/)(^ + 3)] 
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D 



ph 



(5.13) 



An equally good agreement is observed for the function /(p) we consider here, at 
least in the relevant region of p not too large (p < 3), if we use 6 = 1/z/ — 3, and 
Dph = 0.273, Aph = 0.060, see Fig. |. 
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